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The spherically symmetric hedgehog ansatz used in the description of the skyrmion is believed 
to be inadequate for the rotational states such as the nucleon (7 = J = i) and the A (J = J = |) 
due to centrifugal forces. We study here a simple alternative: an oblate spheroidal solution which 
leads to lower masses for these baryons. As one might expect, the shape of the solution is flatter as 
one increases I — J whether the size of the soliton is allowed to change or not. 



I. INTRODUCTION 

When Skyrme first introduced its model a few decades ago |Q] to describe baryons as solitons in a non-linear 
field theory of mesons, the solution proposed was in the spherically symmetric hedgehog ansatz. There are reasons 
to believe that this solution is not adequate for the rotational states such as the nucleon (/ = J = i) and the 
A (/ = J = |) due to centrifugal forces [^[jij. Alternative treatments have been proposed in the past with relative 
success. These approaches generally fall into three classes: (a) The original s^pherical shape of the solution is modified. 
This is usually done by making a global deformation along one or more axisjqj. (b) The size of the skyrmion is allowed 
to change. This analysis also led to the identification of breathing modes [Hpj with excited states of the nucleon and 
A-isobar. Combining deformations (a) and (b), one gets the following scheme: the nucleon's and A- isobar's ground 
states have K = I -I- J = and spherical symmetry. This led to the conclusion that these states are stable against 
quadrupole deformation |^,||. The K ^ states occur in nearly degenerate doublets for all spin: one oblate and 
one prolate, (c) Finally, the shape of the solution itself could be improved by including the (iso-) rotational kinetic 
energy in the minimization of the static Hamiltonian. It turns out that there is no finite static hedgehog solution 
unless one considers massive pions This instability is understood to come from the emission of pions from a 

rapidly rotating skyrmion. Some progress has been made recently to elucidate the connection between skyrmions and 
Feynman diagrams in an effective field theory but it is also interesting to note that the solution were characterized 
by small quadrupole deformations away from the spherical hedgehog ansatz ||9|,p^ . 

In this work, we take a naive approach and propose a simple alternative. Instead of the spherically symmetric 
hedgehog solution, we introduce an oblate spheroidal solution. This leads to lower masses and quadrupole deformations 
for these baryons. Moreover, the shape of the solution is flatter as one increases I = J whether one allows the size of 
the soliton to change or not. 

In the next section, we introduce the oblate spheroidal coordinates. We then proceed as follows: We determine the 
profile / (77) by solving the differential equation obtained from minimization of the static energy. This procedure is 
similar to that of ref. |ll| and one indeed recovers the profile of the hedgehog ansatz in the limit of parameter d ^ 0. 
In Section III, we compute the masses of the nucleon [I = J ~ ^) and of the A-isobar {! = ,]= |). These masses 
get contributions both from the static and rotational energy and will in general depend on the choice of d. The value 
of d for each baryon is fixed by minimizing its mass with respect to d. In the last section we discuss our results with 
respect to scale transformations and check that the above scheme persists. 

II. THE STATIC OBLATE SOLITON 

The oblate spheroidal coordinates (?/, 9, 0) are related to Cartesian coordinates through the expressions 

X = d cosh T] sin 6 cos ip 

y = d cosh r] sin 9 sin (j) (1) 
z = dsmhricos9. 

A surface of constant i] corresponds to a sphere of radius d flattened in the z-direction by a factor of tanh?/. For 77 
small, the shape of the surface is more like that of a pancake of radius d whereas for large 77, one recovers a spherical 
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shell of radius r = Let us note that in the limit d ^ Q, r] ^ oo with r remaining finite, the coordinate system 
becomes 

(i, y, z) = r (sin cos (/), sin 6 sin </>, cos 9) 

which means that it coincides with the spherical coordinate system. The choice of the parameter d determines at 
what scale the "oblateness" becomes important. The element of volume is given by 

dV = -d^ (cosh r/) (cosh^ f] - I + cos^ 9) ■ dryd(cos e)d(t) (2) 

We would like to replace the hedgehog solution for the Skyrme model by an oblate solution. Writing the Lagrangian 
for the Skyrme model 

where Wd^U with U G SU{2). We get the usual expression for static energy density 
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'fJri^UU)-^^Tr{[U,L,f). (4) 
Let us now define a static oblate solution by 

U = e'(^-^)/('') (5) 

where t) is the unit vector fj = . The boundary conditions for the winding number = 1 solution are / (0) = tt 
and / (oo) = 0. Note that this is not a priori a solution of the field equations derived from the Skyrme Lagrangian. 

Using the oblate ansatz U = exp [i [t ■ fj) f [rf]] we get, after a straightforward but tedious calculation and an 
integration over the angular variables 9 and 0, the static energy contributions E2 and E/^ such that 



dVE2 = ^ • - • d?7 (a2i/" + a22 sin' /) 



(6) 



with 

"21 (??) = 2 cosh 77 



"22 (»?) ^ 2 (-2 cosh 77 + (2 cosh' - l) ^ (??)) 

where L(,y)^ In (ffl^) and 

^^^4 - 4^ • ^ • / rf?? ("41/" Sin' / + a42 sin^ /) (7) 
^ 4d Jo 



with 



a4i (??) = 2L (77) 

, , 11 1 , r / XX 2 cosh 77 

"42 (??) = 77 I n~ (2 cosh 77 + L (77)) ' 



2 \ cosh' 77 (cosh' Tj — \) j 

Here, we have expressed the parameters of the model in terms of the constants: 

A=— d^^d. (8) 

It might be required to add a pion mass term to the Skyrme Lagrangian. This term takes the usual form 

2 p2 

^rn = {TrU - 2) (9) 
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leading to the expression 



A 

where e^r = and 

arn {v) = cosh 1] I sinh + ^ 

Minimizing the static energy with respect to / (77) , we need to solve (numerically of course) the non-linear ordinary 
second-order differential equation which reads: 

= 32ehld^ -amsinf + y (-2/"a2i - 2/'a^i + 2a22 sin / cos /) 

+ . ((4 sin3 / cos /) ^42 + (-2/" sin^ / - 2/'^ sin / cos /) ^41 - (2/' sin^ /) a'^^) . (11) 
Ad 

For calculational purposes, we need to set the value of the parameters of the Skyrme Model. and e are first 
chosen to coincide with those of ref. : 

= 129 MeV e = 5.45 = (12) 

= 108 MeV e = 4.84 = 138 MeV (13) 

obtained by fitting for the masses of the nucleon and the A in the hedgehog ansatz. 
The solution near 77 ^ has the form 

/ (77) TT - ai77 

whereas in the limit 77 ^ 00, one recovers the spherical symmetry with, 



fiv) 



2777, 



exp ( (14) 



where ai and k are constants which depend on d and 777,7^. The solutions of differential equation ( |11| ) are presented 
for several values of d and 777^ = in Fig. 0. For smaU d (here d < 0.0001), we get exactly the solution of the 
spherical hedgehog skyrmion. As one increases d, we observe a displacement of the function / (77) and the continuous 
deformation of the soliton from a spherical to an oblate shape. The static energy. Eg — Em^ + E2 + E^, has a minimum 
for d — which is expected since it corresponds to the spherical solution (see Fig. H). 

The masses of the nucleon and of the A-isobar get contributions both from the static and rotational energy and 
will generally depend on the choice of d. We fix the value of d for each baryon by minimizing its mass with respect 
to d. 



III. COLLECTIVE VARIABLES 



Using the oblate solution, we can then compute the masses of the nucleons (/ = J = i) and of the A-isobar 
(/ = J = |). However, several remarks are in order before we go on. When one departs from the spherical symmetry 
of the hedgehog ansatz, it is customary to introduce extra collective variables for isorotation in addition to those 
characterizing spatial rotation since these are no longer equivalent, in general. The spin and isospin contributions to 
the rotational energy are however equal in our case since we use solution (|^) and we are 

only interested in ground states with K = J + 1 = (see Appendix). As a result, we need only consider one set of 
collective variables. 

Let us work in the body-fixed system and assume that the time dependence can be introduced using the usual 
substitution 

U A{t)UA''{t) (15) 
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where A{t) is a time-dependent SU{2) matrix. This transformation leaves the static energy (or mass of the soliton) 
invariant. We can then go on and treat A{t) approximately as quantum mechanical variables. The calculation 
procedure is fairly standard (see ref. for example). 

Using ([l5|) , the Lagrangian gets new terms due to the time dependence of A: 



and 



where 



^4 



dVCi 



32e2 



dV Tr 



Lq, Li 



Lo = AU^A^do (AUA^) 
Following straightforward but lengthy calculations, we get after angular integrations 



(16) 



(17) 



(18) 



and 



where 



Li = -al'Tr 



Li = ^ial'Tr 



■,aU 



iA^A 



Tr 



Tr 



T.A^A 



T.A^A 



(19) 



(20) 



with 
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Aire 
A 
47re 



• 1287r2e4d2 



d 

2^0 
M Jq 



drj cosh rj siv? f Aij 

JO 

drj cosh T] sin^ / [dj sin^ / + Bij /'^j 



An = A22 — ^cosh^ r/ ^ \ ^ (3cosh^ 77 — 4) + coshrysinh^ 77 L (77) 
^33 = 4 cosh'' ^7 — cosh^ V ^ ("^ cosh^ — 3 cosh^ 77 + cosh 77) L (77) 
Bii = B22 = 2 — cosh^ ^ 2 '^'^^^'^ ^ cosh^ 77) L (rj) 

B: 
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2 cosh^ rj + cosh 77 (l — cosh^ r]) L (rj) 



Cii — C22 — 

c, ' 



-1 



(10 cosh 7; - 8 cosh^ ?7 + (4 cosh"* 77-9 cosh^ 77 + 3) i (7;)) 



2 cosh rj 



4 cosh 77 ^ 

((4 cosh^ rj — cosh^ 7; — l) X (r;) — 8 cosh^ 7; + 2 cosh?/) 



(21) 



(22) 



Non-diagonal terms for Aij , Bij and Cij give zero contribution upon (f> integration due to the axial symmetry of the 
solution. 

Let us now consider the quantity 



Tr 



■,A^A 



Tr 



T.A^A 



for an axially symmetric system where a^^ — af^ ^ c?^ and — for i ^ j. We can rewrite 



Tr 



iA^A 



Tr 



r.A^A 



I'J = Tr 



{aA^) + ( 



Tr 



T3 



{a: 



^3 - a") 



(23) 



In terms of the Euler angles 8, $ and ^E", the traces correspond to the expressions 
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(24) 



which finaUy leads to 



Li 



Tr 



{ii + cose$) 



sin"' e$ 



(25) 



+ cos e$ 



(26) 



where 



c = (af + af) 



(27) 
(28) 



Here h and c play the role of principal moment of inertia. 

Quantization of ( p6| ) is straightforward. It indeed represents a symmetrical top with the rotational kinetic energy 
in space and isospace 



E, 














1) 



(29) 



where |J| and |I| are the spin and the isospin respectively and, J3, the z-component of the spin. We have already 
used the relation J3 = —I3 here which follows from axial symmetry. Added to the static energy Eg, it leads to the 
total energy 



J, J 3 



(30) 



Up to now, we have analyzed the rotational and isorotational kinetic energy from the point of view of the body- fixed 
frame. Observables states, however, must be eigenstates of |J|^ , J3, |I|^, /a with eigenvalues J{J + 1), wj, /(/ + 1), 
TO/ where the operators now refer to the laboratory system (as opposed to body- fixed operators in ( p9[ ) and above). 
These eigenstates can be represented by direct products of rotation matrices 



(a,^,7|J,TOj,TO) {p,a,T\I,mj,-m) = D:^^^^ {a,/3,j) dI^^_^^ (p,cr,T) 



(31) 



where (a,/?, 7) and {p,a,T) are, respectively, the Euler angles for the rotation and isorotation from the body-fixed 
frame to the laboratory system. Since we have axial symmetry in the body- fixed system where J3 = —/a, the quantum 
number denoted by m must be opposite in sign in space and isospace rotation matrices. It is convenient to label the 
basis by the sum of the body-fixed spin and isospin, K. The eigenstates of K" are linear combinations of the basis 
states (|3l| ) with Clebsch-Gordan coefficients {J,m;I,~m\K,(S). The explicit calculation of the energy of rotation 
requires in general the diagonalization E^'^^^ . (see ref . |5| for more details) . Since we are only interested in the ground 
states here, i.e. the nucleon and A-isobar, we set K = which simplifies much of the above procedure. 

We proceded with the case TOtt = 0. Numerically, the minimization of the static energy for the spherical symmetric 
ansatz gives Es = ^ ■ [8.20675], = ^ • [8.906] and Ma = ^ • [11.703]. For the oblate spheroidal ansatz, the 

solution for / (77) is found from ( pT[ ) and the parameter d is chosen in order to minimize the mass of the corresponding 
baryon. In general, as d increases, the static energy Eg deviates from its lowest energy value given by the spherical 
hedgehog configuration. On the other hand, oblate configurations have larger moment of inertia which tends to 
decrease the rotational kinetic energy (see Fig. ||). The existence of a non-trivial oblate spheroidal ground state for 
the nucleon and the A-isobar, as it turns out, depends mostly on the relative importance of static and rotational 
energy. 

Our results are summarized in Table |. We find that the ground state for the nucleon is almost spherical but 
nonetheless oblate with d — 0.0013 thus exhibiting a small quadrupole deformation and a slightly lower mass with 
respect to a spherical configuration. For the A-isobar, the oblateness or quadrupole deformation is even more impor- 
tant and accounts for a 4% decrease in mass. We obtain a minimum for the A mass for a value of d = 0.32 with 
Ma = ^ • [11.293]. 
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Since the minimum of the ground state is affected by the oblate shape of the solution, the parameters F^^ and e as 
given in ref. no longer reproduce the quantities they were designed to fit. However, the existence and the form of 
an oblate ground state for baryons depends on the precise value of i^^, e and d through eq. (pi]). Therefore in order 
to fit for Mi^i ~ 939 MeV and A/a = 1232 MeV, we must readjust and e which determine the value of d for the 
nucleon and A-isobar respectively. After several iterations, we find Ft, ~ 118.4 MeV and e = 5.10 with d ~ 0.0014 
(d — 0.40) for the nucleon (A-isobar). 

The numerical calculations for ^ lead to similar conclusions. Starting from input values for F^ and e in ([T^), 
we get a deformation parameter of d = 0.0009 for the nucleon and d — 0.18 for the A-isobar leading to small decreases 
in their respective masses. The deformation parameters here are significantly smaller than what is observed in the 
77i7r = case, which is partly explained by the sensitivity of d with respect to -FV and e. But since the chiral symmetry 
breaking term contributes here (i.e. Em^), this is also connected to the relative importance of the rotational energy 
contribution to the baryon mass and perhaps more importantly to how each contribution depends on d. 



IV. DISCUSSION 



Quadrupole deformations were found previously [p|pC|] in the context of rotationally improved skyrmions. Contrary 
to our variational approach, these solutions involve the minimization of a Hamiltonian which includes the (iso-) 
rotational kinetic energy, i.e. eq. ( pi] ) with contributions from the (iso) rotational kinetic energy. Nonetheless, we 
found that the oblate spheroidal ansatz gives lower energy than the spherical one for baryon ground states. 

Of course, ansatz (^) is not necessarily the lowest energy solution, the latter being obtained in principle by solving 
the integro-differential equation of ref. jl^]. Unfortunately, only large-distance asymptotics of this solution can be 
written in a closed form. Moreover, the most relevant physical quantity here, the mass of the baryons, gets negligible 
contributions from that region and so it is not very sensitive to the exact form of the solution at large distances. Yet, 
it may be interesting to consider how our results compare to large-distance asymptotics of the solution given in ref. 
pO| in which the pion field reads 



^(r,J) 



r— i-oo B 



r 



exp ( — TOTr?") ( J ■ I" ) J 



exp 



— \/ mi — 



J X rxJ 



(32) 



Here j'^ = (Imk)^^ where j'^ and Imk are the Skyrmion classical angular momentum and moment of inertia tensor 
respectively. As one might expect, we recover the hedgehog solution form in the limit — > and r — > oo. For 
7^ 0, the second term in (|3^) dominates which can be interpreted as a swelling of the Skyrmion with the pion 
field pointing in a direction perpendicular to J due to centrifugal forces. Unfortunately, the exact magnitude and 
direction of J can only be obtained by solving the full integro-differential equation. On the other hand, for the oblate 
ansatz the pion field takes the form Tr{rj) = fjf [rj) and coincide with the hedgehog solution in the limit of large 
distances. The magnitude of fcin (|l^) is found by solving ( pi] ) and optimizing for the deformation parameter d. Our 
numerical calculations show that k increases slowly with d, which suggests that the configuration of the Skyrmion at 
large distances is hedgehog-like and swelling for increasing isospin. This is in partial agreement with the qualitative 
features of ( |3^ ) . Be that as it may, we recall that the purpose was mainly to look at possible deformations at middle- 
range distances since this is where energy and baryonic densities are the largest. Ansatz (]^) turned out to be a rather 
simple, intuitive and efficient trial solution. 

It may also be interesting to consider deformations of the oblate skyrmions under scaling of the unitary transfor- 
mations U (r) such that 



U (r) = Uo (pr) 



(33) 



to minimize the total energy of the nucleon and A-isobar. This corresponds to skyrmions which are allowed to 
change in size. Recall that the previous calculations proposed a change in shape (oblate spheroidal vs spherical). The 
treatment is straightforward and indeed very similar to that of ref. [g . In our calculations the scale transformation is 
equivalent to the substitution 



d^i 
p 



(34) 
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Rewriting the expression for the masses in the body-fixed frame as: 



M'"'' = Em^ +E2 + E4+ ^ ' + .,33^ 33, 35 



we see that under the scaling transformation (p3) 



The total energies Mpf{p) and M/^{p), computed in the laboratory system, can be minimized with respect to the 
p parameter, i.e. to the energically favored size of the oblate skyrmion. The results are shown in Table | for both 
the oblate and spherical cases. The baryon ground states are now swelled oblate solutions. Again, one should in 
principle readjust the F^^ and e parameters to fit the masses of the nucleon and A-isobar. It would also be interesting 
to readdress the problem of breathing modes with these oblate skyrmions. This is a problem for further research. 
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APPENDIX: 



The rotational energy for an axially symmetric system is given by 

—D {ojiili + 0)2^2) 

where uj and are the angular velocities in coordinate and isospin space. A, B, C and D are positive quantities cor- 
responding to moments of inertia. The previous expression can be written in terms of body-fixed angular momentum 
J and isospin I: 

= [A (If + II) + B {Jl + J|) + 2D {h,h + /2 J2)] + § 

where li and Ji are the body-fixed components of J and I. Here we have already used the fact that J3 = —I3 because 
of axial symmetry. Clearly, the nucleon and A-isobar ground states are obtained when K = J + I = since the term 
2D (/i Ji -I- I2J2) is then negative. The general expressions for A, B. C and D are rather lengthy and therefore not 
given here (see rcf. However, they become much simpler for a solution of the form (||) in which case A = B = D 

such that the rotational energy of the nucleon and A— isobar ground states are given by 




The spin and isospin contributions to the rotational energy are equal in this case. 
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FIG. 1. Solutions of differential equation ( p] ) for several values of d (F^ = 129 MeV, e = 5.45 and = 0). For d = 0.0001, 
we get exactly the solution of the spherical hedgehog skyrmion. 



FIG. 2. Static and rotational energies for the nucleon as a function of d in units of ^j^. 

TABLE I. Ground states for the nucleon and A-isobar. The results are shown for both the minimum oblate spheroidal 
configuration and the spherically symmetric ansatz for comparison. The values of M"'''^^ are defined according to eq. ( |3C| ) 
whereas M"^"'^(pmin) is minimized with respect to the scaling parameter p (see eq. (js^)). All masses are expressed in units of 
with parameters = 129 MeV, e = 5.45 and = 0. 









Oblate (d > 0) 






Spherical [d — 0) 






d 




Pmiii 






Pmin 




Nucleon 


0.0013 


8.904 


0.868 


8.797 


8.906 


0.867 


8.799 


A 


0.32 


11.312 


0.670 


10.064 


11.703 


0.668 


10.238 
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